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Chirality Selection in Open Flow Systems and in Polymerization 
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As an attempt to understand the homochirality of organic molecules in life, a chemical 
reaction model is proposed where the production of chiral monomers from achiral substrate is 
catalyzed by the polymers of the same enatiomeric type. This system has to be open because 
in a closed system the enhanced production of chiral monomers by enzymes is compensated 
by the associated enhancement in back reaction, and the chiral symmetry is conserved. Open 
flow without cross inhibition is shown to lead to the chirality selection in a general model. 
In polymerization, the influx of substrate from the ambience and the efflux of chiral products 
for purposes other than the catalyst production make the system necessarily open. The chiral 
symmetry is found to be broken if the influx of substrate lies within a finite interval. As the 
efficiency of the enzyme increases, the maximum value of the enantiomeric excess approaches 
unity so that the chirality selection becomes complete. 

KEYWORDS: homochirality, open system, polymerization, enantiomer, rate equation, dynamical chiral 
symmetry breaking 



1. Introduction 

Between two possible stereostructures of organic 
molecules, i.e. a right-handed (D) and a mirror- image 
left-handed (L) form, life on earth has chosen only one 
type: L-amino acids and D-sugars. 1 Various mechanisms 
for the origin of this chiral asymmetry, so called homochi- 
rality, have been proposed, 1 but the predicted asymme- 
try turned out to be very minute, and therefore, it has 
to be amplified to explain the homochirality. 

Frank showed theoretically that an autocatalytic pro- 
duction of chiral molecules with an antagonistic process 
amplifies enantiomeric excess (ee) and brings about the 
homochirality. 2 His theory has had no experimental sup- 
port for a long time. Recently, an example of ee ampli- 
fication has been found in the production of pyrimidyl 
alkanol, 3 and the temporal evolution was explained by 
the second-order autocatalytic reaction. 4 In our previous 
papers, 5-8 we have reported that, in addition to the non- 
linear autotacalysis, recycling of achiral substrates by the 
decomposition of chiral products accomplishes complete 
homochirality in a closed system without any antagonis- 
tic processes. 

The chemical reaction model we proposed 5 involves 
production of a chiral molecule C from achiral molecules 
A and B. With an assumption that there is an ample 
amount of the material B, the reaction is governed solely 
by the concentration a of the substrate A. We adopt the 
convention of R, S representation to denote two enan- 
tiomers of the chiral product C as (ii)-C and (S)-C, 
which arc further abbreviated as R and S, respectively, 
for brevity. The reaction is assumed to include a nonlin- 
ear autocatalytic effect as well as a recycling back reac- 
tion, and the concentrations r and s of chiral species R 
and S develop by the following rate equations; 

r = (ko + kir + k 2 r 2 )a — X r, 
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is = (ko + fcis + k2S 2 )a 



A s, 



(1.1) 



where x represents the time derivative x — dx/dt. The 
reaction system is assumed to be closed so that the con- 
centration a of the substrate A is determined by the 
conservation as a = en — r — s. Here, the constant Cn 
is fixed at the initial time as cq = ao + tq + sq. Reac- 
tion coefficients k\ and k% correspond to those in the 
production processes without autocatalysis, with a lin- 
ear and a quadratic autocatalysis, respectively. Ao is a 
rate coefficient of decomposition process from the prod- 
uct R or S to the substrate A. The major conclusion we 
have drawn out from the model eq.(l.l) is that, only in 
the case with a quadratic autocatalysis and decomposi- 
tion process (k 2 , A > 0), the chiral symmetry breaking 
is possible. 5 

However, the relevance of pyrimidyl alkanol to life is 
suggestive at best, and the problem of the homochirality 
in life is still not yet resolved. One way to address this 
problem is to construct a model that incoporates some 
characteristic features of organic molecules in life. For 
instance, amino acids have the potential to polymerize 
into chain molecules, as peptides and proteins. They act 
as catalysts or enzymes to produce various molecules, 
and among them may be amino acids themselves, al- 
though this expectation has not so far been confirmed 
experimentally. Based upon this hypothesis, Sandars 9 
proposed a "toy model" for the generation of homochi- 
rality such that polymers catalyze production process of 
chiral monomers from achiral substrates. He included 
cross inhibition effect such that the addition of wrong 
handed enantiomer to the polymer halts further polymer- 
ization. 9, 10 The cross inhibition is similar to the mutual 
destruction effect introduced by Frank. In this paper, we 
would like to propose an alternative model for the gener- 
ation of homochirality, exploiting a chracteristic feature 
that reactions take place in an open system; we do not 
require processes of cross inhibition. 

Most of amino acids (19 out of 20 species ) are chiral, 
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but they are composed of achiral molecules as carbon 
dioxides, anmonium, water and so on. These achiral sub- 
strates are supplied from and leave away to the ambience. 
On the other hand, the products, amino acids, polymer- 
ize to form peptides and proteins, which not only func- 
tion as catalysts or enzymes but also are used for many 
other purposes as body construction, metabolism, etc. 
This means that some amino acids leak away from the re- 
action path of catalysts production: The chiral products 
thus flow out from the reaction system in concern. All 
these features lead us to the necessity to study an open 
system. We construct a minimal reaction model includ- 
ing the above ingredients. In a certain limit, we can show 
that open reaction systems are reducible to our previous 
model eq.(l.l) or similar models in a closed system. 5 

In the following section §2, a simple model of chiral 
molecule production under an open flow is shown to be- 
have similar to that in a closed system after a transient 
time. In §3, the polymerization model is studied, and the 
open flow is found necessary to break the chiral symme- 
try for this model. The results are summarized in the last 
section §4. Some details of stability analysis is relegated 
in the appendix. 

2. Chemical Reaction in an Open System under 
Flow 

First we extend our nonlinear autocatalytic system 
proposed previously 5 to an open system under flow, and 
show that a flow plays essentially the same role as the 
recycling process in a closed system. 

As for the origins of the incoming flow, substrate 
molecules A may be produced by some chemical reac- 
tions inside the system, or there may be a supply of A's 
by a flow from ambience. In any case, the incoming flux 
of A is assumed to be constant and is denoted as F. On 
the other hand, dissipation may be due to the outgoing 
flow which takes away part of all the chemical species 
from the region where the reactions are taking place. 
For simplicity the rate of outflow is assumed common 
to all chemical species, and denoted as /io- Combining 
the effects of autocatalytic chemical reaction and in- and 
out-flows, the rate equations are written as 

f = f(r)a - X Q r - fj, r, 

s = f(s)a - X Q s - /i s, 

a = -(/(r) + f(s))a + A (r + s) + F - ^a. (2.1) 

Here f(r) = fco + far + A^r 2 is the production rate of 
R including the effect of autocatalysis, /(s) is the corre- 
sponding one for S, Ao is the rate of back reaction from 
chiral products R or S to the substrate A. It is then 
easy to show that the sum of all the species relax to a 
constant value Co = F/fXo in a relaxation time /1q 1 as 
r + s + a = Co + (ro + sq + ao — c ) exp(— fi t). Thus, af- 
ter the relaxation time Hq 1 , the open system we consider 
reduces effectively to a closed system (1.1) at a fixed to- 
tal concentration Co with a modified back reaction rate 
A = Ao + /io. The difference between the two systems is 
essentially the value of Co; in a closed system it is deter- 
mined as an initial condition, whereas in an open system 
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it is controlled by the in- and out-going flows. There- 
fore, even in the case that the back reaction is absent 
(Ao = 0), the outflow plays essentially the same role as 
the back reaction, and the chirality selection can be re- 
alized in combination with nonlinear autocatalysis. This 
conclusion leads us to propose an experiment of the Soai 
reaction 3 under the flow with an expectation that the 
complete homochirality can be realized. 

The analysis in this section reveals that the generic 
form of the model (1.1) has a wide range of applicability 
to describe the chirality selection in various systems. In 
the following section in the analysis of polymerization 
processes, we encounter a generalized form of this type 
of rate equations 

r = f(r)a — A(r)r 

a = f(s)a - X(s)s, (2.2) 

where f(r) is an effective production rate and A(r) is an 
effective rate of back reaction. The substrate concentra- 
tion a follows another generic equation. Equation (2.2) 
does not contain the cross inhibition term such as — A c rs . 
By suitable conditions for the functions f(r), A(r) and 
other parameters, the chirality selection is realized as is 
demonstrated in the following. Some of the details is de- 
scribed in the Appendix. 

3. Chirality Selection in Polymerization 

We now study the chiral symmetry breaking in the cat- 
alyst polymerization processes. Our model starts from 
chemical reactions of chiral molecule production from 
achiral substrate, as 

A^R, A^S. (3.1) 

Here an achiral substrate A reacts to produce chiral 
species R and S with a rate fco and the reverse reaction 
takes place with a rate Ao- R and S symbolize two enan- 
tiomeric forms of amino acids in L and D forms. We call 
them R and S chiral monomers respectively hereafter. 
Since the production of amino acids requires energy, the 
rate fco is very small. Only in some special environment 
and places, for instance under thunderstorms or close 
to seafloor hydrothermal vent, the energy is supplied to 
increase fco- But if there is only the enhancement of ran- 
dom production by fco, then the racemic mixture of R 
and S monomers results, as is shown in the following. 
Thus, for chirality selection, we have to consider further 
mechanism. 

Chiral monomers are assumed to polymerize if they 
are of the same enatiomcric type as 

ki fej 

R + R% y~ Ri+i j S + Si y~ oj-i-i , (3.2) 

with i = 1, 2, • ■ • . Here fa and \ represent respectively 
the polymerization and decomposition rate of the chi- 
ral polymers Ri+i or Sj+i. Generally, polymers at higher 
levels or macromolecules have some chemical functions. 
When the 7V-mers, Rat and Sat, act as enzymes to re- 
produce the monomer of the same enatiomeric type, the 
catalytic process is described as 

A + R N ^R + R N , A + S N ^S + S N . (3.3) 
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Here k e represents the rate of catalytic production and 
A e that of back reaction. A good enzyme should have a 
large production rate k e compared to kg , that of the non- 
autocatalytic production process. However, the ratio of 
rate coefficients of a forward and a backward reactions 
is specified by the energy difference of initial and final 
states, A and R (or S). Then the back reaction with A e 
should also be enhanced so as to satisfy the following 
relation 

k e k 

For simplicity, further polymerization starting upward 
from the enzyme iV-mers will not be considered in the 
following: namely, the reaction processes (3.2) only up to 
i < N — 1 arc considered here. 

While the whole system is described by the reaction 
processes (3.1-3.3) for a closed system, there are other 
processes for an open system; that is the exchange of 
substrate molecules between the system and the envi- 
ronment E as is illistrated as 

E^A, (3.5) 

Mo 

and the leakage of polymers out of this system into the 
environment as 



rAe. 



(3.6) 



The leakage process does not necessarily mean the out- 
flow of the materials to the ambience but can be con- 
sumptions of polymers for other utilities. 

In the present paper, we consider the simplest case 
such that the dimer is sufficient to present catalytic ef- 
fect, and assume N = 2. Then, the concentrations of 
achiral substrate a and of chiral products up to the 
dimcrs r, r 2 , s, s 2 evolve respectively as 

2fcir 2 



r =koa — Aq?' + k e ar 2 — X e r 2 r 
r 2 =kir 2 - Air 2 - [i 2 r 2 , 
s =koa — Aq.s + k e as 2 — X e s 2 s ~ 2k\s 2 



2Air 2 - /iir, 



\()S i- K e US 2 — /\ e S 2 S — Zfi,'lS -f 2Al,S2 — UlS, 
S 2 =kiS 2 - Ai,S 2 - [l 2 S 2 

a = — 2k a + A (r + s) - k e a(r 2 + s 2 ) + A e (r 2 r + s 2 s) 

+ F-[i a. (3.7) 

These rate equations yield that the total concentration 
of molecules c = a + r + 2r 2 + s + 2s 2 evolves according 
to 



F - fi a - /xi(r + s) - 2/i 2 (r 2 + s 2 ). 



(3.8) 



Let us first consider a closed sytem with F = /Xj = 
0, (0 < i < 2). Then c remains constant c = cq. 
Even with this conservation law, rate equations are still 
too complicated to study the chirality selection analyt- 
ically, since there are many variables coupled nonlin- 
early. Therefore, we utilize a steady-state approximation 8 
such that the dimerization proceeds promptly and one 
may set f 2 = s 2 = 0. Then the steady-state concentra- 
tions of dimers are determined in terms of monomers as 



f-2 = (fci/Ai)r 2 and s 2 = {k\/ '\i)s 2 , and the rate equa- 
tions for the monomers reduce to 



=(fc 



<k 



k e ki 2 A 
——r a - — r) 
Ai k 



k e k\ r 



A, 



0(a 



(3.9) 



with a = cq — (r+s) — 2(r2+s 2 ). Here the relation eq.(3.4) 
is used. Dimer enzymes lead to nonlinear autocatalytic 
effect, but the effect is compensated by the same amount 
of enhanced back reaction. Since fco + {k e ki/\i)r 2 and 
fco + (fc e fci/Ai).s 2 are always positive, the sole possible 
fixed point is racemic: f = s = (fco/Ao)a. Only in the 
absence of the back reaction, Ao = A e =0, the system 
has a fixed line a = 0. There is no chirality selection if 
the system has a fixed line, as was pointed out in the 
previous studies of a closed system. 5,6 

From the above analysis, it is concluded that the poly- 
merization system (our model at least) has to be open 
for the chiral symmetry breaking to occur. To analyse the 
open system theoretically, we again use the steady-state 
approximation by replacing the dimer concentrations by 
their stationary values, r 2 = k\r 2 /(Ai + ^ 2 ), s 2 = 
fciS 2 /(Ai + /i 2 ). The approximation does not alter the 
structure of the fixed points for the system. The reduced 
rate equations for monomers are obtained as 

k e k\ 



<k 



(Ao 



Al + /!2 

A e fci 



r 2 )a 



Al + /i2 

k e k\ 



Mi 



2/^2 k 



(A, 



Ai + fl 2 
A e fci 



s z )a 



Ai + /i 2 



. + M i + T^ S ) S . (3.10) 
Ai + [i 2 Ai + /i 2 

The enantiomeric excess order parameter of the 
monomer is defined by 

T — S 

<h = —r-- (3-11) 



1 \ 
— r)r, 



From cq.(3.10) it is shown to evolve according to 

$1 = A<Pi - B<j)\ (3.12) 



with 



A = B- 



2koa 



B= o k .\ (r + S) \ak e -\ e (r + s)-2^]. (3.13) 
2(Ai + H2) 

One notices that Ao, /io and fi\ do not appear explicitly 
in eq.(3.13), but their effects are implicitly included in 
the values of concentrations, a and r + s. 

As long as the concentrations a and r + s are positive, 
the coefficient B is always larger than A. If A is positive 
asymptotically or at a fixed point, then 4>i approaches 
to a nonzero value as <t>\ = A/B and the chiral state 
will be selected. For example, if A > at racemic fixed 
points (f = s), they are unstable and the system may 
evolve to chiral fixed points. On the contrary, if A < at 
racemic fixed points, they are stable. For //q = ^2 = 0, 
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the racemic fixed point is easily obtained from eq.(3.7) 
and eq.(3.8) as 

f = S -=A ^ = Mi + ^o + (WAi)f 2 _ 
2/Ui' fco + {k e ki/Xi)r 2 

By using the relation eq.(3.4), the coefficient A in 
eq.(3.13) is written as 

A = —k e 

x A e (fc 1 ?^ 2 /A 1 ) 2 - Qi! - 2A )(fc 1 f 2 /A 1 ) + A (A + t ix)/X e 
fc e (fcif 2 /Ai) + k 

(3.15) 

Therefore, A is positive and the racemic fixed point is 
unstable, if the influx F is in the range between the lower 
bound F\ and the upper one F2, given as 



The degree of chiral symmetry breaking is measured 
by the total enantiomeric excess 

4> = Ts - Ss . (3.18) 



F h2 = ^YY^ -2AoTmvT-8Ao/Aii) (3.16) 

with fix > 8Xq. This region where the racemic fixed point 
is unstable is denoted by " chiral" in the F — Ao/fco phase 
space in Fig. 1(a). Parameters are fixed at fco = 0.1, fci = 
1, k e = 10 as well as at Ai = /xi = 1, and the decomposi- 
tion rates back to achiral substrate Ao, A e are varied by 
keeping the relation (3.4) satisfied. At a very large F, the 
racemic fixed point recovers stability, since f = s increase 
and the chirality suppression factor A e (f + s) in B due 
to a finite back reaction rates Ao, e exceeds the chirality 
enhancement factor ak e ~ X e f. Therefore, by decreasing 
the back reaction, the racemic fixed point is expected to 
remain unstable at large F. Actually, at Ao = A e = 0, the 
upper bound F2 diverges and the instability of racemic 
fixed point takes place for all F's larger than the lower 
bound 



^1 = 2/ii 



Aifcp 



(3.17) 



When the back reaction is too strong as Ao > H1/8, 
or when the substrate influx is too strong as F S> i*2, 
there remain plenty of achiral substrates and the ran- 
dom production by fco dominates over the catalytic pro- 
duction with fc e . When the substrate influx is too weak 
as F < F% , there are too few product chiral molecules to 
sustain catalytic chiral enhancement. In those regions, 
racemic fixed point is the only possibility. With a small 
decomposition Ao < /Ui/8 and with a moderate substrate 
influx, states with broken chirality are possible. 

We note that the conclusions of the instability on the 
racemic fixed point remain valid for the full system, irre- 
spective of the approximations. Equipped with these re- 
sults, we carry out numerical integration of the full rate 
equations (3.7) systematically. All concentrations except 
r (or s) are set zero initially, and the initial value of r 
(or s) is varied to obtain the flow trajectories. The flow 
diagram in the phase space of the total concentrations of 
two enantiomers, r E = r + 2r2 and s s = s + 2s2, shows 
that the chiral symmetry breaking actually takes place in 
the region denoted by "chiral" in Fig. 1(a). For example, 
flow trajectories at F = 0.5 with Ao/fco = A e /fc e = 0.6 
in Fig. 2(a) clearly shows the symmetry breaking: Fixed 
points with broken chiral symmetry attract all the tra- 
jectories. 



Within the chiral region, the fixed point value of <p varies 
as the incoming flux F, as shown in Fig. 3(a). Due to 
the symmetry, there is a trivial symmetric branch with 
negative (j> values, which is not shown. Near above the 
lower bound Fx , the ee increases continuously. Here, chi- 
ral fixed points lie close to the racemic one, and the 
parameter A in eq.(3.13) is approximately written as 
A cx (F — Fi)(F 2 — F). Since the F and f are small, 
contributions of dimers r 2 in r y is negligibly small. Thus 

<t> PS 4>1 CX y/F-Fx. 

As F increases, <f> increases as well, and the chiral fixed 
point extends for F larger than F2, where the racemic 
fixed point recovers its stability. In fact, for F between F2 
and F 3 there are two types of stable fixed points, racemic 
and chiral, in the flow diagram, as shown in Fig. 2(b) at 
F = 0.8. If the initial concentrations r and s for the 
numerical integration of eq.(3.7) are close to the racemic 
configuration, the final state remains racemic, since it is 
stable, at least locally, as proven by the stability analysis. 
However, if the integration start from a well developped 
chiral configuration, the chiral fixed point is selected. The 
concentration phase space of r's and s's is thus devided 
into multiple basins of attractions to racemic and chiral 
fixed points. When F is larger than F3, the basins of 
attraction to the chiral fixed points seem to disppear. 

From Fig. 2(b) it is evident that for F between F2 
and F3 there are unstable chiral fixed points (marked by 
white diamonds) between the racemic and stable chiral 
fixed points (marked by black boxes). Therefore, in the 
F — <j) diagram in Fig. 3(a), there is another branch cor- 
responding to these unstable chiral fixed points, denoted 
by a dashed curve: it extends from F3 back to the small 
F direction below the stable chiral branch. The unstable 
branch seems to be connected to the F axis at a point 
F 2 . 

The maximum value of ee in Fig. 3 is rather small. 
This is due to the present choice of a parameter set, 
especially fco, which is chosen rather large on purpose 
so that the fixed points are clearly discernible in figures. 
With a parameter set in Fig. 2(a), for example, the ee 
takes a value = 0.871, but with a smaller value of fco 
as fco = 0.01, the ee increases up to <fr =0.992 and the 
stable chiral fixed points lie almost on two axes. 

Every aspect of the phase boundary and the flow tra- 
jectories remains the same, even if all the out-flow rates 
Hi take finite values, for example [io = fJ-i = /i2 = 1: 
Phase diagram in F versus Ao/fco hi Fig- 1(b) looks simi- 
lar to that in Fig. 1(a). When F takes the value between 
Fx and F2, the racemic fixed point is unstable and the 
chiral symmetry breaks as shown in Fig. 2(c). Between 
F 2 and F 3 , there is a hysteresis such that the station- 
ary state depends on the initial concentrations of r's and 
s's. From the lower bound Fx the ee increases continu- 
ously until it drops suddenly to zero at F3, as shown in 
Fig.3(b), where the expected unstable chiral fixed points 
are omitted. 
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We now consider cases without any direct back reac- 
tion. Xi — for i = 1 to 3. Still, part of chiral products 
is utilized for other purposes and disappears from the 
catalyst production route, such that iVs are finite. By 
fixing /io = fj>2 = 1 and at various /xi, the chiral sym- 
metry breaking is achieved at an incoming flux F larger 
than the lower bound, as shown in Fig. 4(a). One notices 
that there is a minimum in the lower bound of the in- 
coming flux F. With /io = /ii = — 1, the ee increases 
rapidly to the saturation value <p close to unity, as shown 
in Fig. 3(c). With small fco and large k e , the saturation 
to (j) = 1 happens almost instantaneously. 

With finite rates of back reactions, as A^ = fc, for ex- 
ample, the out-going fj,i has to be sufficiently large to 
realize the chirality selection, as shown in Fig. 4(b). The 
symmetry breaking is possible for F between F\ and F3 
with hysteresis above F 2 , similar to the previous cases in 
Fig.f (a) and (b). The ee increases continuously from zero 
at the lower bound Fx, and drops abruptly to vanish at 
the upper bound F3 with hysterisis above F2, as shown 
in Fig.3(d). 

4. Summary and Discussions 

We extend our model of chiral molecule production 
from achiral substrate with nonlinear autocatalysis to an 
open system under flow. Even if the back reaction is not 
explicitly included, the system after a transient period is 
shown to be effectively described by the rate equations 
we derived for a closed system with back reaction, and 
the chirality selection is expected to take place. As an 
ideal system to realize homochirality in laboratory, we 
proposed the Soai reaction under flow. 

Organic molecules in life, amino acids for instance, 
may not undergo nonlinear autocatalytic reaction, but 
they polymerize and acquire some functions. One of them 
is to catalyze various chemical reactions, presumably the 
monomer production process as well. We construct a 
model of chiral monomer production from achiral sub- 
strate with further polymerization to produce enzymes 
which enhance monomer production of the same ena- 
tiomeric type. By assuming for simplicity that dimers 
already act as enzymes, the rate equations of monomer 
concentrations reduce to those with quadratic autocatal- 
ysis in a steady-state approximation. In a closed system, 
however, the enzyme enhances the decomposition process 
as well, and the chiral symmetry is not broken. Only in 
an open system where substrate is supplied and parts 
of produced monomers and polymers are consumed for 
other purposes, the chiral symmetry can be broken if pa- 
rameters take proper values. Influx of substrate should 
be large to yield sufficient chiral products to sustain cat- 
alytic process, but if the back reaction takes place the 
influx should not be too large. Otherwise, an enhanced 
back reaction compensates the chirality enhancement. 

We have considered only a single chiral species for 
simplicity reason, whereas in reality there are 20 types 
of amino acids. If only the same enantiomeric type can 
form stable polymers, we believe that the present mecha- 
nism of chirality selection is still valid for multiple chiral 
species. 

By extending the ladder for polymerization processes 



until the enzyme is produced, the window in parame- 
ter space for the chiral symmetry breaking becomes nar- 
rower. Various disturbances due to spatial and temporal 
variations may tend to narrow this window further. How- 
ever narrow the parameter window may be, it is essential 
that the window is open. After all, it seems to be a ma- 
jority concensus that several hundred million years were 
needed for the birth of life on the earth. 
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Appendix 

A general set of the (reduced) rate equations for the 
concentrations r,s,a is 

r = f(r)a — X(r)r 

s = f(s)a-X(s)s (A- 1) 

c = g(a) — h(r) — h(s) 

c = a + j(r) +j(s) 

where f(r) and A(r) are the effective coefficients of 
the production reaction and the back reaction, respec- 
tively, c is the total concentration and functional forms 
of g(a), h(r), j(r) depend on models but presumably are 
nonnogativc functions. An enantiomeric excess order pa- 
rameter defined as 

V — s 

<t> = — (A-2) 



is shown from eq. (A-l) to evolve according to 

4> = A4>- Bcj) 3 (A-3) 

with 

2/(0)a 



A — B - 



r + s 



(A-4) 



B = \{r + s) 

x f s(f(r) - /(0)) - r(f(s) - f(0)) a X(r) - X(s) -j 
I rs(r — s) r — s J 

The form of the parameter B is so chosen that it remains 
finite for the racemic state r = s as well as at the trivial 
point r = s = 0. For specific form f(r) = fco + fc'iT + fc2^ 2 
and A(r) = Ao + Air + A2r 2 , B becomes simply 

B = ^(r + s){k 2 a-X 1 -X 2 {r + s)}. (A- 5) 

We should not forget that A and B are functions of the 
variable r, s,a and varies as functions of time accodingly. 

For simplicity, let us assume that fsa ^ at a fixed 
point r = f, s = s, a = a, which satisfy 

f(f)a = A(f)f, f(s)a = X(s)s. (A-6) 

Using these equations, one obtains 

B = h? + s)!@* (A-7) 
2 rs 

and A = (f> 2 B. The last relation only shows the consis- 
tency. To analyse the stability of the fixed point, one 
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should first to obtain the fixed point values of a, r, s as 
functions of relevant parameters and finds the conditions 
for the chiral phase A > or for the racemic phase A < 
as is discussed in section 3. 

The other way to analyse the stability is the lin- 
ear stability analysis. As a concrete example, we con- 
sider the stability matrix for the full system eq.(3.7) 
at the racemic fixed point eq.(3.14) supplimcnted with 
f 2 = hf 2 /(Ax + (M 2 ), s 2 = fci.s 2 /(Ai + fi 2 ). Since there 
are five independent variables as (xi, x%, X3, X4, x$) 

■ r 2 , s, s 2 , a) , the linear deviation around a fixed point 
Xi is assumed to be Xi = x\ + a^o exp(rW). The equa- 
tion to determine the eigcnfrcqucncy n is an eigenvalue 
equation |T — Sll\ = 0, where / is the 5x5 dimensional 
unit matrix and the matrix elements of T are given as 
Tij = dii/dxj evaluated at the fixed point. In a matrix 
form, it is expressed as 

det 

(s-n 

2kif 


V 



w 



-k e a + X e r. 



T 








Ai -0 











-s-n 


T 





2kxf 


-Ai - n 


w 


V 


w 



u \ 


u 



-2T- nj 

= (A-8) 
where the nonzero matrix elements are given as 

S — fii + A + 4fcif + A e f 2 

T = 2\i + k e a — A e f 

U = k + k e f 2 (A-9) 

V = A + A e f 2 



This equation is expressed as a product of a quadratic 
and a cubic equations as 

[n 2 + (s + \i)n + SAi - 2kiTr] 
x [-{n + s)(n + x 1 )(n + 2U) 

+ 2fciT7^(f7 + 2U) + 2UV(n + Ai) + 4UWkif] = 0. 

(A-10) 

The instability condition for the racemic fixed point is 
the existence of positive eigenfrequency. As S + Ai > is 
always satisfied, an inequality S\i — 2kiTf < is just the 
condition for the quadratic equation to have a positive 
root and it can be shown that this condition agrees with 
that given in section 3. As for the cubic equation, its 
roots are near — S < 0, — Ai < 0, — 2T < and thus are 
expected to be negative. Numerically we confirmed this 
expectation for reasonable parameter ranges. 
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Fig. 1. Phase boundary between the racemic and the chiral fixed points in the F — \o/ko phase space. Production parameters are fixed 
at fcrj = 0.1, fci = 1, fc e = 10. Other parameters fixed are (a) Mo = M2 = 0, Ai = Ml = 1, and (b) Mo, 1,2 = Ai = 1. 




Fig. 2. Flow diagram in the total concentrations of r and s phase space. Black boxes represent stable fixed points, and white diamonds 
unstable ones. Production parameters are fixed at fco = 0.1, fci = 1, k e = 10. Figures (a) and (b) correspond to Fig. 1(a) with 
decomposition rates Ai = l,Ao/fco = A e /fc e = 0.6 and flow mo = M2 = 0. The in-coming flux of the substrate is (a)F = 0.5 and (b) 
F = 0.8. (c) corresponds to Fig. 1(b) with decomposition rates Ai = 1, Ao/fco = A e /fc e = 0.5 and flow /^o = Ml = A*2 = 1- The in-coming 
flux of the substrate is F = 2. 
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Fig. 3. Fixed point value of the total enantiomeric excess <j> versus influx F. Production parameters are fixed at ko = 0.1, k\ = 1, k e = 
10. Figure (a) corresponds to Fig. 1(a) with decomposition rates Ai = 1, Ao/fco = ^e/k e = 0.6 and out-flow fi± = 1, Mo = = 0. (b) 
corresponds to Fig. 1(b) with decomposition rates Ai = l,Ao/feo = A e /fc e = 0.5 and out-flow fio = fj,i = fi2 = 1. (c) corresponds to 
Fig. 4(a) without decompositions, A^ = and out-flow /iq = fii = fi2 = 1- (d) corresponds to Fig. 4(b) with decompositions Ai = ki and 
out-flow fio = M2 = 1, Ml = 2. 




Fig. 4. Phase boundary between the racemic and the chiral fixed points in the F — phase space. Production parameters are fixed at 
kg = 0.1, k± = 1, k e = 10, and Mo = = !• Rates of back reactions are (a) A^ = 0, and (b) \i/ki = 1 with i = 0, 1, e. 



